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Abstract. In this article combinatorial aspects of normal ordering annihilation and 
creation operators of a multi-mode boson system are discussed. The modes are as- 
sumed to be coupled since otherwise the problem of normal ordering is reduced to the 
corresponding problem of the single-mode case. To describe the normal ordering in 
the multi-mode case for each mode a colour is introduced and coloured contractions 
are considered. A depiction for coloured contractions via coloured linear representa- 
tions is given. In analogy to the single-mode case associated coloured Stirling numbers 
are defined as coefficients appearing in the process of normal ordering powers of the 
number operators. Several properties of these coloured Stirling numbers are discussed. 



PACS numbers: 02.10.Ox 



1. Introduction 



Since the seminal work of Katriel pQ the combinatorial aspects of normal ordering 
arbitrary words in the creation and annihilation operators and b of a single-mode 
boson having the usual commutation relations 

[6,6+] = 66+ -6+6 = 1, [6,6]=0, [&+,&+] =0 (1) 

have been studied intensively and many generalizations (e.g., g-deformed bosons) have 
also been considered, see, e.g., 0[3ISE>l[IE[71[Hl[3[!I2ra 

and the references given therein. For an important discussion and references to the 
earlier literature on normal ordering noncommuting operators see Wilcox [20 1 . For the 
creation and annihilation operator p and / of a single-mode fermion having the usual 
anticommutation relations 

{/,/ t } = // t + / t / = l, {/,/} = 0, {/ t ,/ t } = 

the analogous combinatorial problem does not exist due to the nilpotency of the opera- 
tors, i.e., 

(/t)2 = = / 2 . However, if one considers instead of a single-mode fermion a 
multi-mode fermion (i.e., several sets of operators /j, / ') then interesting combinatorial 
connections to rook numbers exist (note that the general normal ordering problem of a 
single-mode boson also has connections to rook numbers [E]). This was noted by Navon 
[21] even before Katriel demonstrated that normal ordering powers of the bosonic num- 
ber operator, i.e., (b^b) n , involves the Stirling numbers of second kind [lj. To the best 
of our knowledge the normal ordering of creation and annihilation operators of a multi- 
mode boson has not been treated in the literature in a combinatorial fashion analogous 
to the single-mode case, although multi-mode systems are of particular importance in 
applications, e.g., in quantum optics [22J. Clearly, to obtain interesting results the dif- 
ferent modes have to be coupled (i.e., have a nontrivial commutation relation), otherwise 
the problem reduces to the simultanous normal ordering of the individual modes. In the 
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present paper we assume that all modes b\, hi satisfy the commutation relations ([T]) and 
are coupled via [&i,o]] = 1 (all other commutators between operators of different modes 
vanish). In this paper we begin the combinatorial study of the normal ordering of this 
multi-mode boson. It is important to note that the situation is slightly different from 
the one typically considered in quantum optics: There the different modes commute 
and the nontrivial interaction is introduced via the Hamiltonian; here we don't consider 
Hamiltonians but only the coupled creation and annihilation operators. Also, we discuss 
the relations coming from the algebra itself and consider only casually representations 
of the algebra. 

The paper is structured as follows. In Section 2 we discuss the two-mode boson, in 
particular the normal ordering problem. It is shown that the general normal ordering can 
be described with the help of the double dot operation and two-coloured contractions, 
a generalization of the coventional contractions. Some examples are also considered 
and the relation to the normal ordering process of a single-mode boson is discussed. In 
Section 3 we describe a graphical representation of two-coloured contractions using two- 
coloured linear representations. In Section 4 we introduce two-coloured Stirling numbers 
as coefficients appearing in the normal ordering process of the two-mode boson. Several 
properties of these two-coloured Stirling numbers are discussed, e.g., recursion relations 
and some explicit values. In Section 5 we discuss representations of the commutation 
relations of the two-mode boson and give a simple example. In Section 6 it is described 
how the results of the preceding sections can be generalized to the general multi-mode 
boson. Finally, in Section 7 some conclusions are presented. 



2. The two-mode boson 

In this section we discuss the two-mode boson. That means that we consider two sets 
\b\ b} and {a^ ,a} of bosonic creation and annihilation operators which each satisfy (H|) 
and are in addition coupled via the commutation relations 

[a,^] = l, [a,b}=0 (2) 

(together with the adjoint equations [b,a'] = 1 and [a', b>] = 0). Since we will often 
compare this two-mode boson with the case of two noncoupled bosons it is convenient 
to introduce a terminology which will be used in the following. 

Definition 2.1. The coupled two-mode boson (or, briefly, two-mode boson) consists of 
the operators {a, a\b, b^} satisfying the commutation relations 

[a, a] = [a f , a f ] = [6, b] = tf] = [a, b] = [a f , 6 f ] = 0, [a, a)] = [b, 6 f ] = [a, 6 f ] = [b, a)} = 1. 

(3) 

Similarly, the noncoupled two-mode boson consists of the operators {a, a) , b, b^} satisfying 
the commutation relations 

[a, a] = [at, a f ] = [b, b] = [b\ 6+] = [a, b] = [a 1 ", 6 f ] = [a, b+] = [6, a 1 "] = 0, [a, a 1 "] = [6, 6 f ] = 1. 

(4) 

In both cases we speak of the modes a and b (consisting of the operators {a,at} resp. 
{6, and also imagine that the modes have different colours. 

In the following we are interested in properties of the coupled two-mode boson and - 
if not stated otherwise - the commutation relations ([3]) are assumed to hold. 
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Let us introduce the number operators N a and N b as well as the number-difference 
operator by 

N a := a f a, N b := tfb, N d := N a - N b = a) a - b ] b. (5) 

It is straightforward to show that one has the usual relations [N a , a] = —a, [N a , a^] = 
and [N b , b] = -b, [N b , tf] = 6 1 " as well as the "cross relations" [N a , b] = -a, [N a , o 1 "] = a) 
and [A^,a] = —b, [N a ,tf] = b^ . Note that these relations imply that a — b commutes 
with N a and with N b , i.e., [N a ,a - b] = = [N b ,a - b]. Also, one has [N a ,N b ] = 
a)b — b^a. The number-difference operator Nd plays a prominent role in recent studies 
on superconductors, see, e.g., |23j (but note that the two modes discussed there commute 
in contrast to ([2])). 

The normal ordering is a functional representation of two-mode boson operator func- 
tions in which all the creation operators stand to the left of the annihilation operators. 
Although it is irrelevant from a strictly logical point of view (since ab = ba as well 
as aftf = &W) we will assume in the following that in the normally ordered form the 
operators precede the operators and that the operators b precede the operators a. 
Let an arbitrary word F(a, a) , b, tf) be given; a function F(a, , b, tf) can be seen as a 
word on the alphabet {a, a* , b, b^}. We denote by J\f[F(a, a', b, b^)] the normal ordering 
of the function F(a, a', b, b'). Using the commutation relations ([3j) it is clear that its 
normally ordered form J\f[F(a, , b, b^)] = F(a, , 6, tf) can be written as 

Af[F(a,a!,b,tf)]=F(a,a\b,tf)= £ C k ^ n {J) k (tf) l b m a n (6) 

k,l,m,n 

for some coefficients Ckj,m,n an d the main task consists of determining the coefficients 
as explicit as possible. 

Remark 2.2. Let us consider the situation of the noncoupled case ([!]); this is just the 
two-dimensional harmonic oscillator and is the situation considered, e.g., in [23]. Given 
a function F(a, a T , b, ¥), we denote the function which results by deleting all operators 
a, by F a ~ del (b,tf); similarly, we denote the function which results by deleting all 
operators b, b^ by F \a,aJ). Since all operators of the one mode commute with all 
operators of the other mode one has 

F(a, a f , b, &*) = F b - del (a, a))F a - de \b, 6+). 

Clearly, the problem of normal ordering F(a, , b, b^) reduces to the problem of normal 
ordering the single-mode operator functions F b ~ del (a, a)) and F a ~ del (b,b j: ). 

In the case of a single-mode boson one can describe the normal ordering by means of 
contractions and the double dot operation [19J. Let us recall the procedure from [TJ5] for 
a word F(b,b^). A contraction consists of substituting 6 = and b^ = 0^ in the word 
whenever b precedes and deleting all the letters and 0^ in the word. Among all 
possible contractions we also include the null contraction, that is the contraction leaving 
the word as it is. We define C(F(b, b^)) to be the multiset of all contractions of F(b, b^). 
The double dot operation arranges a word it € C(F(b, 6^)) such that all the letters 
precede the letters b (e.g., : b h (b^) 1 : = (¥) b k ). Then one has 

Af[F(b,tf)]=F(b,tf)= £ :tt:. (7) 

7rGC(F(6,fet)) 

We now introduce two-coloured contractions and show that in general the normal 
ordering of an operator function of the two-mode boson ([3]) can be described with the 
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help of two-coloured contractions and the double dot operation. Here the "colours" are 
used to mark the two different modes of operators involved and the use of "colours" will 
become clearer when we consider two-coloured linear representations in Section O 
A two- coloured contraction consists of 

(1) substituting 

(a) a = and a) = 0^ in the word whenever a precedes a', 

(b) 6 = and 6^ = 0^ in the word whenever b precedes 6^, 

(c) a = and = 0^ in the word whenever a precedes b' , 

(d) 6 = and a) = 0^ in the word whenever b precedes o) , 

(2) and deleting all the letters and 0^ in the word. 

Among all possible two-coloured contractions we also include the null contraction, that is 
the contraction leaving the word as it is. We define C^ 2 \F(a, a\ 6, &t)) to be the multiset 
of all two-coloured contractions of F{a,a [ ,b : b^). Following [19] we call a two-coloured 
contraction to be of degree r if precisely r pairs of creation and annihilation operators 
are contracted. 

Example 2.3. Let F(a, at, b, 6+) = aV6 2 6t. Since there are two creation operators 
in the word it is clear that contractions of degree at most two can arise. The null- 
contraction gives trivially a 2 a^6 2 &t. There do exist six contractions of degree one (here 
we have indicated in the brackets which two operators were contracted): 

a6 2 6 t (13), ao) 6 2 (16), a6 2 6 f (23), ao) 6 2 (26), a 2 a f 6(46), a 2 a f 6(56). 

Likewise, there do exist six contractions of degree two: 

6 2 (13, 26), a6(13, 46), a6(13, 56), 6 2 (16, 23), a6(23, 46), a6(23, 56). 

Thus, the set of two-coloured contractions is given by 

C( 2 )(aW6t) 

= {a 2 a)b 2 b J[ , ab 2 b\aa)b 2 , ab 2 b\aa)b 2 , a 2 a)b, a 2 o)b, b 2 , b 2 ,ab, ab, ab, ab}. ' 

The double dot operation arranges a word vr G C (2) (F(a, a 1 ", 6, 6 f )) such that all the let- 
ters a) precede all letters precede all letters b precede all letters a. As an example one 
has : b k a m {b' i ) l {a)) n : = {a)) n {b' i ) l b k a m . It is now possible to describe the generalization 
of (J7]) in the following theorem. 

Theorem 2.4. Let F(a,a^ ,6,6^) be an operator function of the two-mode boson (0). 
The normally ordered form of F{a, a) , 6, 6^) can be described with the help of two- coloured 
contractions and the double dot operation as follows: 

N[F(a,a\b,tf)} = F(a,a\b,tf) = ^ : vr : . (9) 

7reC( 2 )(F(a,at,6,fet)) 

Example 2.5. Let i^a, a^ , 6, 6^) = a 2 a^6 2 6^ from above. The associated set of two- 
coloured contractions is given by ([8]) and it follows from Theorem I2.4I that 

N[a 2 a)b 2 tf] = aWa 2 + 2a f &a 2 + 2a t 6 2 a + 26 f 6 2 a + 26 2 + 46a (10) 

which can also be checked by hand using (J3j) - Let us now come back to Remark 
12.21 and assume that the two modes are not coupled, i.e., satisfy It follows that 
F a - del (b,tf) = 6 2 6t and F b - del {a,a)) = a 2 a) as well as N[a 2 a)b 2 tf] = N[a 2 a))N[b 2 tf) = 
(at a 2 + 2a)(6^6 2 + 26). Multiplying out and sorting the terms yields 

M[a 2 a)b 2 b ] \ = aWa 2 + 2a f &a 2 + 26 f 6 2 a + 46a (11) 
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which clearly differs from the coupled case (jlOp . 

Example 2.6. In this example we again come back to Remark 12. 21 and assume that the 
two modes are not coupled, i.e., satisfy implying that the modes comprise just the 
two-dimensional harmonic oscillator. This implies that [N a , N b ] = which is crucial for 
the following. Let us consider the operator function F(a, a), b, ¥) = N™ a N^ b for some 
powers n a ,n b € N. It follows that N[N™ a N T b lb ] = Af[N^ a ]M[N b nh ] and we may use the 
beautiful result of Katriel PQ that 

Af[K«] = J2S(n a ,k)(J) k a k (12) 

k=0 

(where the S(n a , k) are the Stirling numbers of second kind) to conclude that 

n a n b 

Af[N^K b ] =^Y,S(n a ,k)S(n b ,l)(^) k (b j ) l b l a k . (13) 
fc=0 1=0 

Under these circumstances it is also easy to determine the normally ordered form of 
powers of the number-difference operator since N d = N a - N b and, therefore, N% d = 
Y1U ( n k d )(-^) k K d ' k Nt, implying 

n d n d -k k / \ 

AflK d ] = EEE d )(-l) k S(n d -k,r)S(k,s)(a^(tfyb s ar. (14) 

k=0 r=0 s=0 \ ) 

Given such a two-dimensional harmonic oscillator, a well-known application is to use 
it for a representation of the angular momentum operators. This is the coupled boson 
representation introduced by Schwinger [24] . see also [21 ED] for discussions. Namely, if 

J+ := tfa, J_ := t 6> j z .- i(iv 6 - N a ) = -~N d 

then J_|_ and J_ are the angular momentum raising and lowering operators, respectively 
(and J z is the z-component). 

Example 2.7. Let us consider F(a, cfi, b, tf) = (b r a^a) n = (b r N a ) n . Its normally or- 
dered form can be written as 

AfftVafia)"] = Y,c { : ) (i,j,k)(a^yVa k (15) 

for some coefficients ch\i,j, k). Using a k a^ = a)a k + ka k ~ 1 , one shows 

(a^yVa k b r a^a = (a<y + b j+r a k+1 + (j + r)(a)) i V +r - 1 a k+l + k{<J) i b j+r a k , 
implying the following recursion relation for the Cn \i,j, k): 

ci r ji (h 3, k) =4\i-l,j-r,k-l) + (j + r)c« (i,j — r+l,k — l) + kc^ (i, j -r,k). (16) 

Let us check that this recursion relation reduces for the special case r = to the one 
of the Stirling numbers of second kind (if r = then F(a, cv, b, W) = (a^a) n = iV™ is 
the function of a single-mode boson considered in Example 12. 6p . For r = one has the 
recursion relation 

c ( nUi,j, k) = <4°>(i - 1,3, k - l)+jcW(i,j + l,k-l) + kc^(i,j, k) 

where ch^ (i, j, k) = for j > 0. Thus, defining first T n (i, k) := oh 1 (i, 0, k), this recursion 
relation reduces to T n+ i(i, k) = T n (i — 1, k — 1) + kT n (i, k). From the above definition 
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it is clear that T n (i, k) = if i 7^ k. Thus, only the case i = k remains to be considered. 
Defining S(n, k) := T n (k, k) = oh^ (k, 0, k), the last recursion relation reduces to 

S(n + 1, k) = S(n, k - 1) + kS(n, k) (17) 

which is indeed the recursion relation for the Stirling numbers of second kind. The first 
few explicit values for the coefficients c£ k) are given by 

the only nonvanishing coefficients for n = 1 are 

4 r) (0,r-l,l) =r, cS r) (l,r,l) = 1, 
and the only nonvanishing coefficients for n = 2 are 

4 r) (2,2r,2) = 1, 4 r) (l,2r-l,2) = 3r, 2r, 1) = 1, 

4 r) (0,2r - 1,1) = r, 4 r) (°> 2r - 2 > 2 ) = r ( 2r - !)• 

Given an operator function F(a, a\b, one can associate to it an operator function 
of a single-mode boson by neglecting the colour, i.e., by "putting the sunglasses on". 
To be more precise, we consider on the one hand the unital *-algebra A^ generated 
by the elements {1, a, a', b, b'} having the relations ([3j) (and the involution f) as well 
as the algebra A^ generated by the elements {l,c, c^} having the relations dH). Of 
course, A^ is just the usual Heisenberg-Weyl-algebra. "Putting the sunglasses on" 
means introducing a *-algebra-homomorphism 

SUM : A^^ A^, 

mapping 1 h 1 as well as a\b^ 1— > c' and a, b 1— > c. Given an operator function 
F(a, a), b, b^) € A^ 2 \ we can either first put the sunglasses on and then normal order 
the resulting operator function of a single-mode boson, or we can first normal order 
it and then put the sunglasses on and obtain a normally ordered single-mode operator 
function. It is a beautiful fact that these two results agree, i.e., the order of the two 
operations is irrelevant. 

Theorem 2.8. Let F(a,a) ,b,tf) € A^ be given and let us denote the operation of 
normal ordering in the two algebras by Mi : A® — > A^ for i = 1,2. Then one has 

SUM{M 2 [F(a, a\ b, 6 1 ")]} = Mi[SUM{F{a, a\b, &+)}], 

or, in a more fancy language, 

SUM oM 2 = Mio SUM. (18) 

However, if we want to reduce the combinatorial complexity of normal ordering the 
two-mode boson by mapping it to a single-mode boson and normal ordering there, we 
have to keep track of the original colours somehow. By this we mean that "morally" 
one wants to write M 2 = SUM^ 1 o Mi o SUM (which, of course, makes no sense since 
SUM is not invertible). 

Example 2.9. Let us continue our study of F(a, a) , b, b^) = a 2 a)b 2 b^ € A^ from above. 
"Putting the sunglasses on" yields SUM(a 2 a^b 2 tf) = c 2 Sc 2 S G A^ . A straightfor- 
ward calculation shows that normal ordering this expression gives M\[SUM{a 2 a^b 2 b^)] = 
(cr) 2 c 4 + 6c^c 3 + 6c 2 . On the other hand, we have calculated the normally ordered form 
of a 2 a t 6 2 6 t , i.e., M 2 [a 2 a^b 2 tf] G A^ , in Example [231 with the result CEU]). Applying 
SUM to the right-hand side of CE]) gives SUM(M 2 [a 2 a^b 2 tf]) = (ct) 2 c 4 + 6ct c 3 + 6c 2 
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which agrees with Ni[SUN (a 2 a)b 2 b^)] - as it should according to Theorem 12.81 Note 
that the above-mentioned problem with writing A/2 = SlAN~ l 0A/1 oSUN ', i.e., with the 
"map" SUAf^ 1 , amounts to the problem of reconstructing the right-hand side of (fTOl) 
from (c^) 2 c 4 + 6c^c 3 + 6c 2 - which is not possible without further information. 

Remark 2.10. It is important that in the modes are coupled, i.e., that [a, = 1 = 
[b, gt ] holds true since otherwise Theorem 12.81 - in particular (|18p - would not be true. 
Note that on the right-hand side of (|18p the relation between the two modes a and b 
is irrelevant whereas it is crucial for the left-hand side. Thus, for the equation to hold 
one cannot choose an arbitrary relation between the two modes. As an example, we 
consider again the function F(a, , 6, 6 T ) = a 2 a)b 2 b^ and assume that the two modes are 
not coupled (all operators of one mode commute with all operators of the other mode, 
i.e., dU) holds true). Using QJ, we find that SUN{N 2 [a 2 a)b 2 tf)) = (c t ) 2 c 4 +4c t c 3 +4c 2 + 
Mi[SUN{a 2 aW)]. 

3. TWO-COLOURED LINEAR REPRESENTATIONS 

Let us recall from [19] that contractions in the single-mode case can be depicted with 
diagrams called linear representations. Let us consider a word tt on the alphabet {b, &t} 
of length m, i.e., it = tti • • ■ vr m with itj G {b, b^}. We draw m vertices, say 1, 2, . . . , to, 
on a horizontal line, such that the point j corresponds to the letter ttj. We represent 
each 6 by a white vertex and each letter by a black vertex; a black vertex j can be 
connected by an undirected edge (i,j) to a white vertex i if i < j (but there may also 
be black vertices having no edge). Importantly, the edges are drawn in the plane above 
the points. This is the linear representation of a contraction. 

Now, we generalize this to the coupled two-mode case ([3]) by introducing two-coloured 
linear representations. We assume that two colours c a and c& have been chosen which are 
associated to the modes a and b. Let us consider a word 7r on the alphabet {a, a> , b, b^} of 
length to, i.e., tt = 7Ti ■ • • 7r m with ttj E {a, a T , 6, ¥}. We draw m vertices, say 1, 2, . . . , to, 
on a horizontal line, such that the point j corresponds to the letter iXj. We represent 

(1) each letter a by an empty circle of colour c a , 

(2) each letter a) by a filled circle of colour c a , 

(3) each letter b by an empty circle of colour q>, 

(4) each letter b^ by a filled circle of colour q,. 

A filled circle j can be connected by an undirected edge (i,j) to an empty circle i if 
i < j - independent of the colours of the circles. As above, the edges are drawn in the 
plane above the points. This is the two-coloured linear representation of a two-coloured 
contraction. 

Example 3.1. The two-coloured linear representations for the function F(a, o) , b, b') = 
a 2 o)b 2 b^ from above are shown in Figure [H Note that we have - due to technical reasons 
- represented the different colours by different symbols: The mode a is represented by 
empty of filled black boxes, whereas the mode b is represented by empty or filled circles. 

4. TWO-COLOURED STIRLING NUMBERS 

In this section we introduce two- coloured Stirling numbers by generalizing the obser- 
vation (|12p of Katriel pQ . The most natural analogue to (|12p is to consider the normally 
ordered form of either N£ a N£ b or (N a N b ) n . 
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□ □ BOO* 

1 2 3 4 5 6 

□ rr"*"""i o o**"» 

1 2 3 4 5 6 

123 nr^ 

□ tT^. o <j^% 

1 2 3 4 5 6 

Figure 1. The two-coloured linear representations of the two-coloured 
contractions of the word aaa^bbb^ . 

Definition 4.1. Let n,n a ,n b be natural numbers and let the operators {a,a^,b, b^} 
satisfy the commutation relations §3$) of the coupled two-mode boson. Then the two- 
coloured Stirling numbers of first type S nat n b (i, j, k,l) are defined by 

M[N n aN n b] = J- S na ^ 3 ,k,l){^) l (b^b k a l (19) 
i,j,k,l 

and the two-coloured Stirling numbers of second type S n (i,j,k,l) are defined by 

M[(N a N b ) n ] = J2 S n (iJ,k,l)(^y(b^b k a l . (20) 
i,j,k,l 

Remark 4.2. Let us come back to Example 12.61 Assuming that the modes are not 
coupled, i.e. holds true, we have determined J\T[N™ a N^ b ] in (|13|) . implying that in 
the noncoupled case 

Sn a ,n b {i,j,k,l) = S(n a ,i)S(n b ,j)5 it i5j, k . 

Let us start with the two-coloured Stirling numbers of first type. It is clear that these 
numbers reduce to the conventional Stirling numbers of second kind if n a = or n b = 0. 
More precisely, if n b = then j = = k and i = I and we have the identification 
Sn a ,o(h 0> 0, i) = S(n a , i). Similarly, if n a = then i = = I and j = k and we have the 
identification 5o, ni) (0, j, j, 0) = S(n b ,j). From the definition it is clear that due to the 
two-coloured contractions the number of creation and annihilation numbers are equal 
for all summands, i.e., only summands with i+j = k + l are not vanishing. An expansion 
shows that N2 a N£ b = (a^) na (tf) nt >b nb a na + • • • where the dots mean lower order terms. 
Thus, S na:nb (n a ,n b ,n b ,n a ) = 1 and S na! n b (i,j,k,l) = if i,l > n a or j,k > n b . Let us 
consider some explicit examples for small exponents, beginning with n a = 1 = n b . One 
has 

N a N b = a)b ] ba + o)b, 
so that the only nonvanishing Stirling numbers S\ t i(i,j, k, I) are 

51,1(1,1,1,1) = 1, 5i,i(l, 0,1,0) = 1. 

Similarly, one has 

N%N b = (o + ) 2 6 t 6o a + 2{a)) 2 ba + cJtfba + a)b, 
so that the only nonvanishing Stirling numbers S2 t i(i,j, k, I) are 

52,1(2,1,1,2) = 1, 52,i(2,0,l,l) =2, 5 2 ,i(l, 1,1,1) = 1, 52,1(1,0,1,0) = 1. 
Since 

N a N^ = a\b ] fb 2 a + 2aW + a)b ] ba + a)b 



a □ ^b o o • tf" □ b o o -» 

123456 123456 



□ rf^*^B O O • 

1 2 3 4 5 6 



□ □ ■ tf"""" O ^"""^ □ □ ■ O 

123456 123456 
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the only nonvanishing Stirling numbers £l 2(2, j, k, I) are 

£1,2(1,2,2,1) = 1, £1,2(1,1,2,0) =2, 5-1,2(1,1,1,1) = 1, £i,2(l,0,l,0) = 1. 
Our final explicit example is given by 

N 2 N 2 = (at) 2 (6t)2 6 2 a 2 + 4(a t ) 2 6 t 6 2 a + (at) 2 6 t 6a 2 + 2(a t ) 2 6 2 

+2(a t ) 2 6a + a){tf) 2 b 2 a + aW + 2a t 6 t 6a + a*b, 

showing that the only nonvanishing Stirling numbers £2,2 (*> j,k,l) are 

£ 2 , 2 (2,2,2,2) = 1, £ 2j2 (2,l,2,l) =4, £2,2(2,1,1,2) = 1, £2,2(2,0,2,0) =2, 
£ 2 ,2(2,0,1,1) =2, £2,2(1,2,2,1) = 1, £2,2(1,1,2,0) = 1, £2,2(1,1,1,1) =2, 
£ 2 ,2(1,0,1,0) = 1. 

Using 

a ] a{a)) l {b ] yb k a l = (a)) l+1 {tfyb k a l+1 + j{a)) i+x {tfy- x b k a l + i{a)) l {tfyb k a l , 
one finds the following recursion relation 

Sn a +l,n b (i,j,k,l) = S na!nb (i-lJ,kJ-l)+(J+l)S na>nb (i-l,j+l,k,l)+iS na!nb (i,j,k,l). 
Similarly, using 

(a)j{tfyb k a l b% = (a ] ) l {tfy +1 b k+1 a l + k{a)) i {b^yb k a l + Z(a t )X6 t )^ fc+1 a'- 1 , 
one finds the analogous recursion relation 

Sn a ,n b +1 (i,j,k,l) = S na!nb (i,j-l,k-l,l)+kS na> n b {i,j,k,l) + (l+l)S na>nb (i,j,k-l,l+l). 

Considering the first recursion relation for = (with j = = k, i = I and identifying 
Sn a ,o(,h 0, 0, i) = S(n a ,i), see above), this reduces to £(n a + l,i) = S(n a ,i — l) + iS(n a ,i) 
which is the recursion relation (|17p of the Stirling numbers of second kind. In the same 
fashion the second recursion relation reduces for n a = also to the recursion relation of 
the Stirling numbers of second kind. 

Let us now consider the two-coloured Stirling numbers of second type, S n (i,j,k,l). 
From the definition it is clear that £o(0, 0,0,0) = 1 and that all other numbers with 
n = vanish. Considering n = 1 shows that Si(i,j,k,l) = Si t %(i,j,k,l). Thus, the first 
interesting case is n = 2. Here one has 

(N a N b f = (at) 2 (6t)W + 3(at) 2 6t 6 2 a + (a t ) 2 6 t 6a 2 + (at) 2 6 2 + (a t ) 2 6a 

+a t (6 t )Va + a^tffba 2 + aW + Aa^tfba + a)b. 

Note that (N a Nb) 2 ^ N 2 N 2 . Thus, the only nonvanishing Stirling numbers £2^, j, k, I) 
are 

£2(2,2,2,2) = 1, £2(2,1,2,1) =3, £2(2,1,1,2) = 1, £2(2,0,2,0) = 1, 
£2(2,0,1,1) = 1, £2(1,2,2,1) = 1, £2(1,2,1,2) = 1, £2(1,1,2,0) = 1, 
£2(1,1,1,1) =4, £2(1,0,1,0) = 1. 

From their definition it is also clear that S n (i,j, k, I) vanishes if i + j 7^ k + I. We now 
summarize the properties of the two-coloured Stirling numbers in the following theorem. 

Theorem 4.3. The two-coloured Stirling numbers of first type S nat n b (i,j,k,l) satisfy 
£o,o(^i> k, I) = 5i,o^j,o^fc,o^/,o an d the following recursion relation with respect to n a 

(i, j, k, I) = S na)rib (i — l,j, k,l—l) + (j + l)S na>ni (i - 1, j + 1, k, l)+iS na , nb (i, j, k, I) 

and nb 

(i,j,k,l) = S najnb (i,j-l,k-l,l)+kS na>nb (i,j,k,l) + {l+l)S na)nb (i,j,k-l,l+l). 
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Furthermore, one has S natnb (n a ,nb,rn,,n a ) = 1 and S na)n , b (i, j,k,l) = unless i,l < 
n a,j,k < nf, and i + j = k + I. The two-coloured Stirling numbers of second type 
S n (hj,k,l) satisfy So(i,j,k,l) = 5j,o^j,o^fc,o^/,o an d ^ e recursion relation 

S n+l {i,j,k,l) = S n (i - 1, j - 1, k - 1, 1 - 1) + (j + l)S n (i - l,j,k - 1,1) 

+jS n (i - k,l-l) + (j + l) 2 S n (i + 1, k, I) + iS n (i, j -l,k-l,l) 
+iS n (i,j — l,k,l — l) + i(2j + l)S n (i,j, k, I) + i(i + l)S n (i + 1, k, I). 

Furthermore, one has S n {n,n,n,n) = 1 and S n (i, j, k,l) = unless i,j,k,l < n and 
i + j = k + I. 

Proof. The properties of the S nat n b (i,j,k,l) were shown above. The only nontrivial 
property which remains to be shown is the recursion relation of the S n (i,j,k,l). A 
straightforward but tedious computation shows that 

a i atfb(a i ) i (tfyb k a l = (a t )* +1 (6 + ) J ' +1 6 fc+1 o ,+1 + (j + l){a)) i+1 {tfyb k+l a l + j (a f ) i+1 (6 1 J b k a l+1 

+3 2 {a ] ) l+1 {b ] y- l b k a l + i{a)y{b ] y +l b k+l a l + i{a)) l {tfy +1 b k a l+1 
+i(2j + l)(cJy(tfyb k a l + i(i - l){a}y- x {tfy +l b k a l . 

This yields directly the asserted recusion relation for the S n (i,j, k,l). □ 

The recurrence relations in the statement of Theorem I4.3I can be written by partial 
differential equations as follows. Let us define (as we said I = i + j — k) 



T na ,n b (u,V,w) := S na , nb (i,j,k,i + j - k) 



u l v j w k , 



u l v j w k . 



T n (u,v,w) := ^2 S n (i,j,k,i + j - k) 

Multiplying the recurrence relations in the statement of Theorem I4.3I by u l v^w k and 
summing over all possible i, j, k > we get that 

T na+ i >nb (u,v,w) =u (l + ^ + T natrib (u,v,w), 

T na)nb+l {u,V,w) = W [V + U-§^ + V-§^ + (1 - W)^ T najnb (u,V,w), 

T n+1 (u,v,w) = u 



w{v + i) + {vw + v + i) (e. + a.) + v ( a_ + ^ , 



9\ 2 



T n (u,v,w), 



for all n a ,nb > 1 and n > with the initial conditions T^ifu, v, w) = uw(v + 1) and 
Tq(u, v, w) = 1. The above recurrence relation for T n+ i(u, v, w) can be written as follows. 
Let 

T = T(x; u, v, w) := 2_. T n (u, v, w)x n . 

n>0 

Multiplying the recurrence relation for T n (u,v,w) by x n and summing over all n > 
yields the rather complicated looking partial differential equation 

, Nm , fdT dT\ (d 2 T d 2 T d 2 T\ 

T = 1 + XUW[V + ljT + XU[VW + V + 1) — h — + XMl) -^-j + 2 + 



9u ch; / I <9u 2 Such; <9t; 2 y 

Recall that the conventional Bell numbers may be defined as a sum over the Stir- 
ling numbers of second kind, i.e., B(n) = J2k^( n ^k), and are given explicitly by 
5(0) = 1,5(1) = 1,5(2) = 2,5(3) = 5,5(4) = 15,... After having introduced and 
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discussed some properties of the two-coloured Stirling numbers we may also introduce 
the associated two-coloured Bell numbers of first resp. second type by 

B^(n a ,n b ) := £ S na>nb (i,j,k,l), B®(n) := £ S n (i,j,k,l). 

i,j,k,l i,j,k,l 

From the explicit values for the two-coloured Stirling numbers given above we find the 
explicit values 

£(2) (0,0) = 1, 5 (2) (1,0) = 1, 5 (2) (0,1) = 1, 5 (2) (1,1)=2, 
£ (2) (2, 1) = 5, 5 {2) (1,2) = 5, 5 {2) (2,2) = 15. 

Similarly, one has 

£ (2) (0) = 1, B {2 \l) = 2, B {2 \2) = 15. 

Thus, a comparison with the above-mentioned values of the conventional Bell numbers 
suggests that B^ 2 \n a ,nb) = B{n a + n b ) as well B^ 2 \n) = B(2n). This is indeed true: 

Theorem 4.4. The two-coloured Bell numbers B^ 2 \n a ,nb) of first type and B^ 2 \n) of 
second type are given by the conventional Bell numbers B(n) as follows: 

B^(n a ,n b ) = B(n a + n b ), B^(n) = B(2n). (21) 

Proof. Let us show this for the two-coloured Bell numbers of first type; the proof is 
the same for those of the second type. From the definition of the two-coloured Stirling 
numbers of first type one has M 2 [N^N b n "] = Ei,j,k,i S na ,n b (i, 3, k,l){a^y(tfyb k a l where 
we have indicated with the subsript the algebra in which we normal order. We now want 
to make use of Theorem 12.81 Let us first apply SUM to the above equation, resulting in 

SUM{M 2 [N^N^}} = ]T S na , nb (i,j,k,i+j - fc)( c t)*+ic^ 

i,3,k 

where we have denoted by c, the annihilation and creation operator of the single-mode 
case and have also used that k + l = i + j. Now, we use (]18p and write the left-hand side 
of this equation as Mi[SUM{N^ a N^}] = Mi[N^+ n t>]. Using (HI, we have obtained 

"£S(n a + n b ,l)(J) l c l = ]T S na>nb (i, j, k,i +j - k)(J) i+ ic l+ i. 
I i,j,k 

Setting c = 1 = yields the desired equation B(n a + n b ) = B^ 2 \n a , n b ). □ 

Remark 4.5. In this particular case the operation of "putting the sunglasses on" (i.e., 
mapping to the single-mode case via SUM) proved to be successful - in contrast to the 
case of the two-coloured Stirling numbers. Clearly, the reason for this is that the subtle 
distribution of summands in the two-mode case - which gets lost upon maping to the 
single-mode case - does not matter here since a sum over all possible terms is taken. 

5. A REPRESENTATION OF THE COMMUTATION RELATIONS 

Let us recall that in the single-mode case one has the following representation of the 
canonical commutation relations ([I]): a) i— > x,a i— > ■4- (acting as multiplication resp. 
differentiation operator on an appropriate set of functions). Introducing for the second 
mode the corresponding representation b^ i— > y,b i— ► ^ one obtains a representation of 
the noncoupled two-mode boson dH). 
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We now construct a representation of the coupled two- mode boson (|3|). For this we 
make an ansatz 

d d 
a = q{x,y) — , a f = s(x,y), b = p(x,y) — , ¥ = t(x,y) (22) 
ox dy 

where p, q, s, t are in the moment arbitrary (smooth) functions and the operators act 
on a space D of (smooth) functions f(x, y) which is also to be determined. Inserting 
this ansatz into the commutation relations ([3]) yields conditions on the functions of the 
ansatz. For example, the relation [o, b] = yields the condition 

= [a, b]f = abf - baf = q^ - p± , 

and, therefore, 

dq df dp df 
^ dy dx ^ dx dy 

(where we have used = q^)- In the same fashion one considers the remaining 
commutation relations of ([3]). As a result, one has the following set of five equations for 
the functions p, q, s, t and /: 

ds dt dt ds dqdf dpdf 

ox oy ox oy oy ox ox oy 

To solve this set of conditions, we now make some further assumptions. As a first step, 
we assume that s = t. Then three independent conditions remain from (I23h . namely 

ds ds dq df dp df 

= 1 ' Pa~ = 1 ' Pa~a~ = 1a~^~- 

ox dy dy dx ox oy 

If we make in addition the assumption that p = q, these equations reduce to 

q— = l q— = l q^L^L = q^L^L (24) 
dx ' dy ' dy dx dx dy 

Thus, ff = f|- The structure of the equations in (|24|) shows the following: If we choose 
the functions / in such a fashion that = and if we furthermore choose q such 
that || = || then all equations are satisfied provided that g|| = 1. Consider the 
functions q(x,y) := e x+y and s(x,y) := e _ ( a; +^). They are smooth and satisfy all the 
conditions. If we now restrict the functions / to be smooth and satisfying = 
then all conditions of (|24p and, therefore, of (|23]) are satisfied. Thus, we have found a 
particular representation of the commutation relations ([3]). Let us summarize the above 
observations in the following theorem. 

Theorem 5.1. Let D := {/ G C°°(R 2 ) | = §£} and let 

a = e *+y° i Q t = e -(*+i/) h = e x+v±L h \ = e -{x+v). (25) 
dx dy 

Then the operators a,a<,b,b* act in D and satisfy the commutation relations (0j of the 
coupled two-mode boson. 

Proof. Formally, we have checked by the above construction that choosing these opera- 
tors satisfy the commutation relations ([3|). The choice of D assures us that interchanging 
the derivatives of / in the above calculation was legitimate. The final point to check is 
that the operators act in D, i.e., that Of G D for / £ D and O E {a, a', b, b'}. This is 
straightforward. □ 
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Remark 5.2. With hindsight, the above result is not too astonishing. Note that the 
restriction ^ = ^ on the functions / £ D is extremely severe, reducing the above 

example essentially to the single-mode case since a) = 6' and af = bf for all / E D 
due to the restriction! Thus, it would be interesting to find other, more complicated 
representations of the commutation relations ([3]). 

Let us consider the representation given in Theorem 15 . 1 1 fur ther . The number operator 
N a is given by iV a = o)a and, therefore, represented by since N a f = e~( x+v * > e x+v J^f = 
The same argument shows that Nb is represented by J^. However, since the two 
derivatives coincide in D, one may conclude for this representation that N a = Nf,. 
Stretching the physical interpretation to its limit, we see that this is some kind of "di- 
agonal" representation where the numbers of the two modes always coincide. It follows, 
in particular, that [N a , Nf,] = so that all conclusions derived in Example 12.61 hold true. 
Recall that we discussed in Section 2 that one has directly from the commutation rela- 
tions (i.e., independent of the representation) [N a , Nj,] = a^b — b^a. Inserting (|25p yields 
[N a ,N b ] 

— dx dy w -hich indeed vanishes on every / £ D. This is a nice consistency 

check. 

6. The general multi-mode boson 

In this section we generalize the preceding sections to the general case of a multi-mode 
boson. For the following we assume that n modes are involved, i.e., there are n pairs of 
creation and annihilation operators a\,cn (1 < i < n) each satisfying (pQ) and where the 
cross commutation relations are given for all 1 < i < j < n by [aj,aj] = 1, [ai,dj] = 
(together with the adjoint equations [oj,a|] = 1 and [o|,at] = 0). 

Definition 6.1. The coupled n-mode boson is given by the annihilation and creation 
operators {a\,a\, . . . , a n , a^} satisfying for all 1 < i, j < n the commutation relations 

[a t , aj ] = 0, [at, oj] = 0, [a*, a]] = 1. (26) 

A general operator function F(a%,a\, . . . ,a n ,a^) for this multi-mode boson may be 
interpreted as a word on the alphabet {ai, a\, . . . , a n , a^}. The normal ordering is a 
functional representation of multi-mode boson operator functions in which all the cre- 
ation operators a] stand to the left of the annihilation operators a{. Although it is 
irrelevant from a strictly logical point of view (since OjOj = ajUi as well as ajaj = aja|) 
we will assume in the following that in the normally ordered form one always has that a\ 
precedes ol whenever i < j and that the operators Oj precede the operators aj whenever 
i > j. Thus, a normally ordered form consists of a sum of terms of the form 

(4) h (4r • • • (4)* , «# -a?af. 

We have chosen this particular form of arrangement since it is very convenient for taking 
adjoints (the adjoint of such a term has again this form). As above, we denote by 
Af[F(ai, . . . , ajj] the normal ordering of the function F(a\, . . . , aJJ. 

In analogy to above we now introduce n-coloured contractions: An n-coloured con- 
traction of the word F(a\, • • • , a}J consists of 

(1) substituting = and aj = 0^ in the word whenever aj precedes aj for any 
l<i,j<n, 

(2) and deleting all the letters and 0^ in the word. 
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Among all possible n-coloured contractions we also include the null contraction, that 
is the contraction leaving the word as it is. We define C^ n '(F(ai, - ■ ■ , a)j) to be the 
multiset of all n-coloured contractions of F(a\, . . . , ajj. Again, we call an n-coloured 
contraction to be of degree r if precisely r pairs of creation and annihilation operators 
are contracted. 

The double dot operation arranges a word it € C^ n \F(ai, • • • , ajj) such that 

(1) the letters a\ precede the letters aj whenever i < j, 

(2) the letters aj precede the letters aj whenever i > j, 

(3) all letters a\ precede all letters aj. 

It is now possible to generalize Theorem 12.41 as follows. 

Theorem 6.2. Let F(a±, . . . ,a n ) be an operator function of the coupled n-mode boson 
IW\). The normally ordered form of F(a\, . . . , a n ) can be described with the help of 
n-coloured contractions and the double dot operation as follows: 

Af[F( ai ,...,4)} = F( ai ,...,4)= Yl ( 27 ) 

7reC(«)(F(oi,...,al)) 

In analogy to the two-mode boson where we introduced the algebra A^ we may 
here introduce the corresponding unital *-algebra A^ generated by the operators 
{l,ai,...,aj l } satisfying the above commutation relations (j26|) . Clearly, the operation 
of normal ordering M = M n is a map M n : A^ — > A^ n ' ■ As in the two- mode case we 
have here a *-algebra-homomorphism SUM : A^ — > A^ mapping aj i— > c and a\ i— > 
for 1 < i < n; we also call it "putting the sunglasses on". We now have the following 
generalization of Theorem 12.81 

Theorem 6.3. Let n be an arbitrary integer. Consider the algebras A^ defined above 
and let us denote the operation of normal ordering in these algebras by J\[ n ■ A^ — > A^ . 
Then the operations of normal ordering and "putting the sunglasses on" commute, i.e., 

SUM o M n = Mi o SUM. (28) 

We may also introduce in close analogy to Section H] the general n-coloured Stirling 
numbers of both types as follows. 

Definition 6.4. Let (m, . . . , n n ) £ N n , m £ N and assume that the operators {aj, aj} 
for 1 < i < n satisfy the commutation relations (|26p . The n-coloured Stirling numbers 
of first type S n , ,n n {^1 j • • • > ^Tij 3ni • • • 

,ji) are defined by 

^[K!---K:}= E S ni ,..., n Ji 1 ,...,i n ,j n ,...,^)(alr i ---(aJ i r«a^---af (29) 

il,...,i n 
31, ■■■tin 

and the n-coloured Stirling numbers of second type S n (i\, . . . ,i n ,jn, ■ ■ ■ , Ji ) are defined 
by 

M[(N ai ---N an r}= £ S m (i 1 ,...,i n J n ,...J 1 )(a\) h ---(ai) i «ak---a?. (30) 

il,...,i n 

These n-coloured Stirling numbers can then be discussed like the ones of the two- 
coloured case (explicit values, recursion relations, etc.). Note that we could also de- 
fine less symmetric Stirling numbers which result by "clustering" some of the modes. 
From this point of view the n-coloured Stirling numbers of first type do not cluster the 
modes, i.e., each mode has its own exponent whereas the n-coloured Stirling numbers 
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of second type are the opposite case since there is one exponent for all modes. In be- 
tween are the aforementioned less symmetric versions. For example, one could consider 
N[(N ai ■ ■ ■ N af ,) ll (N ak+1 ■ • • N an ) 1 ' 2 ] and define associated Stirling numbers as above. In 
analogy to the two-mode case one can also introduce the associated n- coloured Bell 
numbers of first resp. second type by 

-B (n) (ni, . . . ,n n ) := ^ S , „ lj ... ) „„(n, . . . ,i n ,jn, ■ ■ ■ ,h), 

il,...,%n 

B^ n \m) := ^2 S m (ii,...,i n ,j n ,...Ji). 

i±,...,i n 

31 in 

Using Theorem 16.31 one may show the analogue of Theorem 14.41 

B (n) (n l7 ...,n n )= B(m + ■■■ + n n ), B^(m) = B(nm). 

Let us now introduce the general n- coloured linear representations to depict n-coloured 
contractions, thereby generalizing the construction of Section El We assume that n 
colours Ci associated to mode a\ for 1 < i < n have been chosen. Let us consider a word 
7T on the alphabet {a\, . . . , a^} of length m, i.e., tt = tt\ ■ ■ ■ ir m with ttj G {a\, . . . , a^}. We 
draw m vertices, say 1,2, ... ,m, on a horizontal line, such that the point j corresponds 
to the letter ttj. We represent 

(1) each letter ai by an empty circle of colour Cj, 

(2) each letter a\ by a filled circle of colour q. 

A filled circle j can be connected by an undirected edge to an empty circle i if 
i < j - independent of the colours of the circles. As above, the edges are drawn in the 
plane above the points. This is the n-coloured linear representation of an n-coloured 
contraction. 

A generalization of the representation discussed in Theorem 15. II to the case of n modes 
is easy to find: 

Theorem 6.5. Let D n := {/ G C°°(R n ) | §§- = §§- for 1 < i, j < n} and define the 
operators ai, a\ for 1 < i < n by 

n — r Xl+-+X„ ® t _ p -{xi + -+x n ) 

a % -e ^ a,-e 

Then the operators ai , a\ for 1 < i < n act in D n and satisfy the commutation relations 
[2b]) of the coupled n-mode boson. 

7. Conclusions 

In this paper we began discussing the combinatorial problem of normal ordering words 
in the annihilation and creation operators of the coupled two- mode boson. It was shown 
that this problem cannot be reduced in a straightforward way to the well-known case 
of a single-mode boson. On the more conceptual side we have introduced two-coloured 
contractions and have shown how the normally ordered form of an arbitrary word in the 
creation and annihilation operators can be expressed as a sum over the two-coloured 
contractions, generalizing the corresponding result for the single-mode case. We also 
described the two-coloured linear representations as a nice depiction of two-coloured 
contractions. In the single-mode case one has the beautiful observation of Katriel that 
normal ordering powers of the number operator involves the Stirling numbers of second 
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kind. In view of this we were led to define in the two-mode case the two-coloured 
Stirling numbers as coefficients appearing when normal ordering powers of the number 
operators of the two-mode case. These two-coloured Stirling numbers are a natural 
generalization of the conventional Stirling numbers and share several properties with 
them. We derived several of these properties, in particular recursion relations and some 
explicit values. Clearly, a more detailed study of these numbers and their combinatorial 
interpretation is desirable. The associated two-coloured Bell numbers were introduced 
and it was shown that they can be expressed through the conventional Bell numbers. 
All of the above aspects of the two-mode boson can be generalized to the case of n 
modes. We sketched most of these points in the last section where we introduced n- 
coloured contractions, n-coloured linear representations and n-coloured Stirling numbers 
and showed how the results mentioned above for the two-mode case can be generalized to 
the n-mode case. Apart from the points already mentioned it would clearly be interesting 
to find nontrivial representations for the commutation relations of the coupled n-mode 
boson, in particular Fock-space like representations. 

Acknowledgements. The authors would like to thank Simone Severini for useful dis- 
cussions. 
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